Nematic order in square lattice frustrated ferromagnets 
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We present a new scenario for the breakdown of ferromagnetic order in a two-dimensional quantum 
magnet with competing ferromagnetic and antiferromagnetic interactions. In this, dynamical effects 
lead to the formation of two-magnon bound states, which undergo Bose-Einstein condensation, 
giving rise to bond-centered nematic order. This scenario is explored in some detail for an extended 
Heisenberg model on a square lattice. In particular, we present numerical evidence confirming 
the existence of a state with d-wave nematic correlations but no long range magnetic order, lying 
between the saturated ferromagnetic and collinear antiferromagnetic phases of the ferromagnetic 
J1-J2 model. We argue by continuity of spectra that this phase is also present in a model with 
4-spin cyclic exchange. 

PACS numbers: 75.40.Cx 75.10.-b, 75.10.Jm, 



The search for a true quantum "spin liquid" — a quan- 
tum magnet which remains disordered at the very lowest 
temperatures — in dimension greater than one has been 
central to research on quantum magnets for more than 
three decades. Following Anderson 1|, most models of 
quantum spin liquids proposed to date have been based 
on frustrated antiferromagnetic (AF) interactions. The 
resulting spin liquid states involve strong singlet bonds 
between spins, which give rise to a gap in the spin exci- 
tation spectrum (for a review, see e.g. 0). 

None the less, the best characterized experimental re- 
alization of a quantum spin liquid is believed to occur 
in two-dimensional films of solid 3 He, where the inter- 
actions between spins are predominantly ferromagnetic 
(FM), and the resulting state is gapless 0, 0]. This 
raises the interesting question of whether the breakdown 
of long ranged FM order offers a new route to a spin 
liquid ground state ? 

The purpose of this paper is to demonstrate that a 
gapless spin liquid can indeed occur in a simple two- 
dimensional model of quantum spins on a square lattice 
with predominantly FM interactions. To this end we 
present a set of analytic and numerical results which show 
how a new type of "bond-nematic" state emerges from 
the ruins of FM order in an extended Heisenberg model 
on a square lattice. Our work is inspired by a number 
of new FM square lattice compounds [2, and has close 
parallels in earlier work on spin chains |jfl. The special 
case of the multiple spin exchange model on a triangular 
lattice — which is believed to describe the magnetism 
of two-dimensional solid 3 He — will be dealt with else- 
where ||- 

While these results are rather general, for concreteness 
in this Letter we consider a square lattice 5=1/2 frus- 
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FIG. 1: Phase diagram for the spin-1/2 J1-J2-K Heisenberg 
model on a square lattice, in the octant where Ji < is FM 
and both J2 > and K > are AF. Solid lines in the J\—J% 
and J\-K planes denote known Neel phases. A single nematic 
phase interpolates between the FM occurring for large FM Ji 
and these two, different, AF ordered states. 

trated ferromagnet defined by the Hamiltonian 
H = 2J\ Sj • Sj + 2 Ji Si ■ Sj 

+K £ ( p i234 + ^1234) - SI (1) 

(1234) i 

where counts nearest neighbor bonds, (ij}2 counts 
next-nearest neighbor bonds, and P1234 performs a cyclic 
exchange of spins on a single square plaquette. We re- 
strict ourselves to the case where J\ < is FM and the 
other two interactions J2 > and K > are AF. 

The essence of our results, summarized in Fig. ^ and 
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FIG. 2: Schematic phase diagram of the FM J^-J^—K ex- 
tended Heisenberg model in applied field, showing how the 
two-magnon instability crosses over from 1 st to 2 nd order m 
applied field and is eventually succeeded by a conventional 
one-magnon instability. The inset shows representative spec- 
tra for I st and 2 nd order transitions into d-wave nematic state 
(high spin states only). 



Fig. [21 is as follows — the competition between FM and 
AF interactions leads to a large "accidental" degeneracy 
in the one-magnon spectrum of the FM phase, which 
then becomes highly susceptible to the formation of two- 
magnon bound states. These condense, and give rise to 
a new singlet phase lying between the classical FM and 
known Neel phases of the model. This phase has a well 
defined Goldstone mode, but no long-range magnetic or- 
der in the conventional sense. We introduce an appropri- 
ate order parameter for this phase, and identify it as a 
d-wave n-type "bond-nematic" state, distinct from con- 
ventional "n-type" and "p-type" nematic states known 
to occur in other closely related models 0, ^J. From 
numerical exact diagonalization calculations we are able 
both to confirm that this phase exists, and to establish 
that it interpolates between the two well-studied limits 
of Eq. Q, the square lattice J1-J2 model for K — > and 
the square lattice J\-K multiple spin exchange model for 

J 2 -^o. 

We now turn to the analysis. For large FM Ji the 
ground state of Eq. is a uniform saturated FM. In 
the FM J\-J% model, for K = this give s way to 2- 
sublattice collinear Neel AF at large J2 In the FM 
J\—K model, for J 2 = 0, the corresponding Neel state is 
an orthogonal 4-sublattice structure 0]. Our analysis 



begins with a very simple observation — at a classical 
level, the transition between FM and Neel states is 2 nd 
order and takes place where \ J\ \ = 2{K + J2). However, 
exactly at this transition, the one — magnon spectrum of 
the FM phase posses entire lines of zero modes [l3j . This 
"accidental" degeneracy has profound consequences for 
magnetism. 

Firstly, it leads to large zero point fluctuations which 
shift the first instability of the FM to somewhat weaker 
AF coupling. We have performed numerical exact diag- 
onalizatiqn (ED) of Eq. for clusters of 16, 20, 32 and 
36 spins [lj| . The evolution of the ground state energy 
indicates that a I s * order transition from saturated FM 
to singlet ground states occurs for | Ji| s» 2.5(K + J2). 

Secondly, these fluctuations also destabilize the neigh- 
boring Neel phase. Semi-classical estimates for the the 
square lattice J\-J% model suggest that the sublattice 
magnetization of its collinear Neel phase vanishes for 
|Ji| ~ 1.97J2 0| • These calculations almost certainly 
underestimate quantum effects. 

These arguments provide good a priori reason to be- 
lieve that the extended Heisenberg model Eq. sup- 
ports a new spin liquid phase for a finite range of pa- 
rameters between its FM and Neel phases. However they 
give little insight into what this phase might be. 

In order to better answer this question, we introduce 
magnetic field h and examine the nature of the first insta- 
bility of the fully saturated paramagnet as the magnetic 
field is reduced. The situation is summarized in Fig. [21 
For parameters well within the classical Neel phases of 
Eq. , the first instability of the saturated state occurs 
at h = B c i and is controlled by the lowest one-magnon 
excitation in applied field. This instability is against a 
conventional two-sublattice canted AF state. 

In the proposed spin liquid region, however, the one- 
magnon spectrum is highly degenerate. In the extreme 
case, J% = 0, K = |Ji|/2, the one-magnon dispersion 
vanishes altogether, so there is no net energy gain in flip- 
ping one spin at any q and the single flipped spins are 
always localized. On the other hand, pairs of flipped 
spins on neighboring sites can propagate coherently un- 
der the action of the cyclic exchange operator P1234, and 
so gain kinetic energy. It is therefore reasonable to ask 
whether the transition out of the FM phase controlled by 
such two-magnon bound states ? 

Both the one-body problem of a single flipped spin- 
1/2 in a saturated FM background |FM) and the two- 
body problem of two interacting flipped spins can be 
solved exactly. We have calculated the energy of both 
states in the proposed spin liquid region and find that, 
at a critical value of field h = B C 2, two-magnon bound 
states of the form \<j>) = J2ij ^ijS^ S~ |FM) with d-wave 
symmetry are gapless and give the first instability, while 
one-magnon states have higher positive energy — see 
Fig. [21 The binding energy of the two-magnon state is 
of purely kinetic origin. In zero field, the two-magnon 
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bound states first become stable in the J\—J% model for 
J 2 /|Ji| = 0.408 (c.f. [3), and in the J x -K model for 
K/\Ji\ — 0.364. These values should be compared with 
the numerically determined boundary of the FM phase 
at (J 2 + K)/\J\\ ss 0.4. In the special case J2 = 0, 
K = 1^1 1/2) the two-magnon bound state has the com- 
pact form 4>ij = -^[Si t j- ei - Sij- e2 ], where ei = (1,0) 
and e2 = (0, 1) connect neighboring sites on a square 
lattice. 

What then happens below B c2 ? In the case of the FM 
J\—K model for J2 = 0, the answer is quite simple. Two- 
magnon pairs undergo Bose-Einstein condensation. For 
this set of parameters there is a weak, repulsive interac- 
tion between two-magnon pairs, so the transition is 2 nd 
order. This is evident from ED spectra — the locus of 
the lowest lying state in the high-spin sector has a con- 
vex curvature as a function of S, so the magnetization of 
the system must evolve continuously at B C 2 — see the 
inset to Fig. [3 Even spin sector states only appear in 
the lowest states in magnetic field. This period two is 
the evidence of condensation of two-magnon pairs. 

In the FM J1-J2 model for K = the situation is 
a little more complicated. The net interaction between 
two-magnon pairs is attractive, and the transition is I s * 
order. Again this is evident in ED from the concave cur- 
vature of the high spin states, which signals a jump in 
magnetization at saturation field. A similar 1st order 
transition appears in the J±-K model for J2 = well 
close to the FM phase boundary. However in all of these 
cases the structure of the low lying states is essentially 
the same — even spin sector states always have lower 
energy than odd spin sector states, and show an alterna- 
tion between the s-wave A\ and d-wave B\ irreps of the 
square lattice space group C± v . We return these points 
below. 

Where bound pairs of magnons condense they give rise 
to a new form of bond-nematic order. This is to say 
that, while the orientation of individual spins remains 
undetermined, the traceless rank two tensor 

O^in^) = \{S?S? + SfSf) - ^(Si ■ Sj) (2) 

does exhibit long range order |l6| . The matrix elements 
of C Q/3 (r%,rj) are connected to the magnon pairing op- 
erator through the relation S~S~ = O xx - O yy - 2iO xy . 
Physically, one can think of this tensor operator as reveal- 
ing the order "hidden" in the spin-1/2 quantum spin- 
liquid by projecting into a symmetrized spin— 1 Hilbcrt 
space of bond variables with long range correlations. 

The order parameter Eq. J2J) is distinct from the site- 
centred n-type nematic order seen in certain models with 
S = 1 0, and from the p-type chiral order found in the 
multiple spin exchange model with AF interactions |l0j . 
However a very similar order parameter was previously 
introduced in the context of one-dimensional frustrated 
spin chains [3. 
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FIG. 3: Spectrum of extended Heisenberg model for N = 36 
spin cluster with Ji = —1, J2 = 0.4, K = 0, showing the 
marked alternation in energies of the lowest lying states in 
even and odd spin sectors. The inset shows the double peak 
structure of the heat capacity cy(T) for clusters of N = 16 
and N — 20 spins and the same set of parameters. 



Having established the unconventional nature of the 
state near saturation, we now turn to the nature of the 
ground state in the absence of applied field. Proof that 
this is a spin liquid follows from the finite-size scaling 
properties of ED spectra. In Fig. [21 we present the spec- 
trum of a cluster of N = 36 spins for J\ = — 1, J 2 = 0.4 
and K = 0, very close to the boundary with the satu- 
rated FM phase. The signature of long range Neel order 
in such a spectrum would be the existence of a set of 
"quasi-degenerate joint states" (QDJS) which form the 
N = 00 ground-sate, presents in every spin sector (up to 
S ~ \/N), with energies scaling as -Eqdjs ~ jj^ and 
are well separated (at finite N) from the lowest excita- 
tions of the continuum (magnons with energies scaling 
as ~ The QDJS specific to (zr, 0) collinear AF or- 

der comprise one A\ and one B\ state for even S, and 
a (twofold degenerate) E level for S odd [lg. The gaps 
to the low lying states in each spin sector should evolve 

as A ~ 4f\oi 7=1, with coefficients a and 3 determined 

by the spin stiffness and spin wave velocity of the Neel 
state, and thus diplay a negative curvature when plotted 
against 1/N [iH]. 

It is immediately clear from the spectrum in Fig. |3 
that the energies of the odd spin states, which include 
the E symmetry states, are well separated from those of 
the even spin states. Thus the spectrum is not compat- 
ible with the expected Neel order. The unconventional 
nature of the ground state is also reflected in the finite 
size scaling of gaps, shown in Fig. The gaps in the 
even spin sectors scale to zero in the expected manner. 
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FIG. 4: Finite size scaling of the energy gaps of extended 
Heisenberg model with Ji = —1, J2 = 0.4, K = 0, for Ai and 
Bi symmetry states (lower) and E ones (upper) in the even 
and odd spin sectors, respectively. 



But the gaps of E symmetry states in odd spin sectors 
do not — in fact they appear to scale to a finite value 
A odd ~ 0.4. 

This is conclusive proof that the FM J1-J2 model does 
not exhibit Neel order in the immediate neighborhood of 
its FM phase. But we can draw a still stronger conclu- 
sion. The spin even-odd oscillation found in high spin 
states persists down to low spin states, and the A\ and 
B\ states found in even spin sectors are precisely the 
same states which emerge from the two magnon instabil- 
ity of the saturated FM in applied magnetic field. The 
existence of this set of QDJS is the signature of nematic 
order with d-wave symmetry. Thus the nematic order 
found near to saturation persists down to zero magnetic 
field. 

Such nematic order has a degenerate ground-state 
composed of even spin states. It has gapless Goldstone 
modes analogous to the lowest magnons states of the 
Neel state, present in the 5 = 1 sector, scaling as h=. 
Results for a prototypical nematic correlation function 
E^^ei^^r + e.)) (a = 1,2) in the ne- 
matic ground state are presented in Fig. El Nematic cor- 
relations exhibit a "stripey" character, and exist in the 
whole system. Because of the underlying d-wave sym- 
metry, correlations on parallel bond have positive sign, 
while correlations on perpendicular bonds are negative. 

The finite temperature properties of the FM J1-J2 
model provides further insight into the bond-nematic 
state. In an inset to Fig. [3] we show the heat capacity 
calculated for the same exchange parameters. It exhibits 
two clear energy scales — a broad structure extending 
to T ~ \Ji\ which reflects the formation of S = 1 ob- 
jects from the paramagnetic "soup" , and sharp peak at 
T ~ 0.1|Ji|, where these S = 1 objects start to order 
nematically. 




FIG. 5: Bond-nematic correlations of an extended Heisenberg 
model with Jj = -1, J 2 = 0.4, K = 0, for a cluster of N = 32 
spins, showing the d-wave character of the nematic phase. 
Correlations are measured relative to the reference bond 1-2. 



Although it is hard to put a precise boundary on the 
stability of the (it, 0) collinear AF state, our ED stud- 
ies suggest that Neel order begins to break down in the 
FM Jt-Ja model for J 2 - 0.6-0.7| J x \. We also have per- 
formed ED calculations for parameter sets interpolating 
between the FM J x -J 2 model for K = and the FM 
multiple spin exchange model for J 2 = 0. These confirm 
that the nature of the state bordering the FM phase does 
not change. Thus a single nematic phase interpolates be- 
tween the two different Neel order parameters in these 
two limits, and the FM phase for large |Ji|, as shown in 
Fig. ID 

To conclude — in the frustrated square lattice FM's 
considered in this letter, the saturated FM ground state 
undergoes a first order transition into a gapless spin- 
liquid state with bond-nematic order. In applied mag- 
netic field, the transition becomes second order and can 
be related directly to the condensation of two-magnon 
bound states. 
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